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Abstract 



We analyse the chiral symmetry in the random ±J XY model on a N x 2 square lattice with 
periodic boundary conditions in the transverse direction. This "tube" lattice may be seen as a two- 
dimensional lattice of which one dimension has been compactified. In the Villain formulation the 
discrete- valued chiralities or charges associated with the plaquettes of the lattice decouple from the 
continuous degrees of freedom. The difficulty of the problem lies in the fact that the chiralities 
interact through the long range "strong" one-dimensional Coulomb potential - which increases 
linearly with distance - as well as through an exponentially decaying "weak" interaction. By 
comparing the ground state energies for periodic, antiperiodic, and reflecting boundary conditions 
in the longitudinal direction, we show that the chiralities and the XY spins have the same zero-T 
correlation length exponent, whose exact value v c = 0.5564... we determine. The equality of 
these correlation lengths even in the presence of long range chirality-chirality interactions lends 
support to the view that chiral glass order cannot be sustained without simultaneous spin glass 
order. 



1 Introduction 



Wc present a study of the interplay between the spin variables and the chiral variables (chiralities) in 
the ± J XY spin glass. The former correspond to the continuous rotational symmetry of this model, 
and the latter to its discrete chiral symmetry (i.e., the invariance of the model Hamiltonian under 
reflection of all the spins with respect to a reference axis), first pointed out by Villain g, ||. Below 
the lower critical dimension, dz, which is believed to be greater than 2 led, the correlation lengths 
associated with the chiralities and with the spin variables diverge as T~ Ua and T~ v % respectively, at 
the zero-temperature (T) critical point. The question of the relation between the two types of variables 
has become of interest following speculations by Kawamura and Tanemura j^] , by Ray and Moore || 
and by Kawamura [Q, prompted by Monte Carlo simulations, that below de, the two correlation 
lengths are different, with v c > v s . This suggests that the chiralities will order more easily than the 
spins in higher dimensions. Consequently, above di there would be a regime of dimensions with long 
range chiral glass order, but without conventional spin glass order. This possibility receives intuitive 
support from the idea that discrete symmetry leads to long range order more easily than continuous 
symmetry does. 

Two recent publications H ^| address this issue analytically. Both these studies, just like the 
Monte Carlo work of jjj, consider the finite size scaling of the ground state energy differences 
between periodic (P), antiperiodic (AP), and reflecting (R) boundary conditions. In one of them, 
Ney-Nifle and Hilhorst || transform the two-dimensional XY ± J spin glass on a finite N x M square 
lattice into a grand-canonical Coulomb gas problem of which, as is well-known, the logarithmically 
interacting charges represent the chiral variables. The charges must take half-integer values on the 
frustrated plaquettes and therefore cannot vanish even in the ground state. In the case of uncorrelated 
disorder, the plaquettes are randomly and independently frustrated with probability ^, and it is not 
possible to find the ground state explicitly. For that reason, the subsequent analytic treatment of 
remains restricted to the example of a rectangular array of frustrated plaquettes with randomly 
distributed intercolumn distances. In this example, the authors find no evidence for a chiral correlation 
length diverging faster than the spin correlation length. By a heuristic argument they extend this 
conclusion to the case of uncorrelated ± J disorder. 

In an earlier investigation, Ney-Nifle et al NHM considered the random ± J XY model on a one- 
dimensional ladder lattice, again in the Coulomb gas representation. This problem is exactly solvable, 
or nearly so, for general disorder, and the conclusions drawn from it are fully coherent with those from 
the two-dimensional model ||. However, this model suffers from the drawback that, in the Coulomb 
gas language, it has only exponentially decaying electrostatic interactions (for reasons explained in 
that work) , so that one may wonder if an essential ingredient of the difficult two-dimensional problem 
has not been lost. 

In the present work, we reconcile the requirements of exact solvability and truly long-range inter- 
actions between the chiralities by studying the ± J XY spin glass on a N x 2 lattice which is periodic 
both in the longitudinal and the transverse direction. We work again in the Coulomb gas represen- 
tation, and apply different boundary conditions. In section 2, we show that, on this two-dimensional 
lattice with one compactified dimension, the electrostatic interaction decomposes into two compo- 
nents. The first one is a "strong" or charge-charge interaction; it is nothing but the one-dimensional 
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Coulomb potential, which increases linearly with distance. The second one is a "weak" interaction: it 
acts between transversely oriented "dipoles" and decays exponentially with distance. We shall refer to 
them as the Coulomb and the dipolar interaction, respectively. The Coulomb gas representation of the 
XY model Hamiltonian involves, in addition to these two interactions, two supplementary "global" 
terms that couple the system's total electric dipole moment to the boundary conditions imposed on 
the Hamiltonian. These extra terms have drawn a certain attention in the recent literature jlo[ [|. [ttj. 
and they play again an important role here. 

We are not able to solve the ground state problem for the complete Hamiltonian. However, we are 
able to conclude that in the large N limit, whatever the boundary conditions, the ground state is one 
of the infinitely many ground states of the "strong" Coulomb interaction combined with one of the 
global terms. The details of the proof (of largely technical nature) of this fact are given in Appendix A. 
This set of ground states consists, roughly speaking, of charge configurations in which the long-range 
Coulomb interaction is screened away as much as possible by the formation of longitudinally oriented 
dipoles, as exhibited in section 3. The degeneracy within this set is lifted by the weak interaction 
and by the second global term, which are therefore responsible for the selection of the ground state of 
the full Hamiltonian and for the energy differences between P, AP, and R boundary conditions. Even 
though we remain unable to say which member of this set is selected as the true ground state, we are 
able to describe (in section 4) the domain walls and domain wall energies involved in passing from one 
boundary condition to another. Using the relation between the correlation lengths and the finite size 
scaling exponent of the ground state energy differences, we conclude in the final section 5, for the first 
time within an XY spin glass with random ± J disorder and having a nontrivial long-range interaction 
between its chiralities, that the spin and the chiral correlation lengths diverge, for T — ► 0, with the 
same exponent v. We determine its exact value, v = log7 1 (3 + 2\/2) = 0.5564. . ., in Appendix B. 
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2 The tube: a compact ified two-dimensional lattice 

In this section we shall exhibit the Hamiltonians of a random XY model on a N x 2 square lattice with 
periodic boundary conditions (PBC) in the transverse direction and successively periodic, antipcriodic, 
and reflecting boundary conditions in the other direction. We call this lattice a tube (see figure I). It 
can be viewed as a two-dimensional lattice of which one dimension has been compactified. 

First we shall recall the same model on the more general N x M lattice and then specialise to the 
tube. The effect of the compactification on the interaction will thereby become clear. 

We consider a random ± J XY model where the spins are two-component unit vectors whose angles 
fa (with a reference axis) take values in (— vr, tt]. Two nearest-neighbour spins, fa and <j>j, have an 
interaction energy — Jcos(fa — <j)j — itij), where J is a constant, and the 7Ty are quenched random 
variables that take the values 




, with probability j 
7T , with probability 5 
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The partition function is 
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Figure 1: Chiralities on a tube lattice. The bonds on the dashed line are identical with the ones on the lower solid line. 
The bonds have quenched disorder variables it and the plaquette centres frustration variables p associated with them. 
According to the definition (2.4), one has, for example, P( XQ ^\) = ('"'(xq,!) + ^o+i — n ( x o,2) ~ 7T x^)/C^ 7r )- 



The sum in the exponential in (2.2) runs over all nearest-neighbour bonds of the periodic lattice 
with the convention that in < i,j > the site j is to the right of i (for a horizontal bond) or above 
i (for a vertical bond). In our notation the site vectors i = (i Xl i y ) have half-integer components 



1 3 

2 ' 2 ' ■ 



2N-1 



and L 



1 3 

2 ' 2 ' ' 
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Since we are interested in the ground state properties of the model, we shall replace ( |2.2[ ) by the 
corresponding Villain expression, which is believed to be equivalent to (2J2) in the large-/? limit jl], ||), 
and is easier to analyse. The Villain partition function is 

Y[d<Pi e <s,j> » ( 2 - 3 ) 

"* i {«„} 

where the ny are additional dynamical variables. These riy are integers and the sum on them ensures 
that the integrand has period 2n in <pi — <f>j. In the following, we set J = 2. 

For each plaquette of the lattice, we define a frustration variable p r , with r — (x, y) a vector with 
integer components x — 1, . . . ,N and y — 1, . . . , M that labels the centres of the plaquettes, 



Pr 



<i,j> 



2n 



(2.4) 



where the sum is restricted to the bonds that define the plaquette r. In ([Q|), = — 1 or 1 depending 
on whether one runs through the corners of the triangle (ijr) clockwise or counterclockwise. The 
frustration variable is integer for nonfrustrated plaquettes and half-integer otherwise. 

In (2/3) one can integrate on the continuous degrees of freedom. The algebra (cf. U ||, §]) includes 



the transformation from the variables to the new discrete variables q r called the "chiralities" of the 
plaquettes. The chirality q r runs through all integers (half-integers) when p r is integer (half- integer). 
One shows that the chiralities interact via a Coulomb interaction (which is why they are also called 
"charges" ) and that they satisfy the neutrality condition 







(2.5) 
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Recently, Ney-Nifle and Hilhorst Q (see also (8) ) extended the mapping of the XY Hamiltonian 
onto a Coulomb gas Hamiltonian by including all the finite size corrections on a N x M lattice with 
various boundary conditions. We will now adapt their results to the tube lattice, for which a simplified 
notation is defined in figure 1. 



2.1 Periodic boundary conditions 

We shall first consider the N x 2 system with periodic boundary conditions [PBC] in the longitudinal 
direction. We denote its partition function by Z r . Starting from the more general model Z v ||, see 
equation (^^), we change variables from mj to the chiralities q r which allows to perform the Gaussian 
integration on the first set of variables, fa. Including all numerical prefactors in Zq, one gets fl 



{q r } n,m \ r / 



(2.6) 



where <$(•, •) denotes the Kronecker delta. The additional dynamical variables n and m run over all 
integers and the q r take integer or half-integer values, as mentioned above. The Hamiltonian 7i P , 
which will be the starting point of our considerations, reads explicitly |j] 

, / N N \ 2 

+ ^ 2 N (m + i £ x(q (Xtl) + q (x , 2) ) + ^ 

+ k 2 J2 qrq r >U N:2 (r - r') . 

t,t' 

We will briefly discuss its meaning. The first two terms are due to the finite system size. They 
represent a coupling of the horizontal and the vertical component of the total electric dipole moment, 
respectively, to the quenched disorder. In the third term, Un.m(R) is the interaction between two 
charges 

1 » e i(Xk X +Yky) _ 1 



with R = (X, Y), k x = 0, . . . , 27r ^ 1 ^ and k y — 0, . . . , 2?r ^ 1 ^ . The asterisk indicates that 
the term (k Xl k y ) = (0, 0) is left out of the summation. 

In d = 2, Un,m (N,M — > oo) is the two-dimensional Coulomb interaction which varies as a 
logarithm at large distances jl|, ^]. For the tube, we will see in what follows that the compactification 
leads to a decomposition of Un,2 m to two parts: a one-dimensional Coulomb interaction that increases 
linearly with distance and an exponentially decreasing interaction, which is a remnant of a two- 
dimensional dipole-dipole interaction. The appearence of the linear Coulomb interaction and its 
competition with the dipolar interaction makes the model interesting. 

To separate these two interactions in Un,2, we combine the two chiralities of a column x as 

it = Q(x,i) + q(x,2) , , 2 ^ 

qx = q(x,i) - q{x,2) ■ 
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Introducing g+ and q x in (2/7) and evaluating (2J3) for N — > oo in these new variables, we get 



7i P = 



A 



^(x,l) 

2tt 



A r 



AT 

71-2 E QxQx' U p( x - x ') 

x,x' — l 
N 

+ E <£<ii>u?(x - x') . 

x,x' — l 

We find that the charges g+ interact via the long-range periodised Coulomb potential 



+ 2vr 2 iV m + £ E a;<£ 



J_ V 

2Af 



2AT ^ s i n 2f*£1 



N -> oo, fixed, £ if < 1 



(2.10) 



(2.11) 



If |X| is negligible with respect to TV, Up(X) is the usual one-dimensional Coulomb interaction, linear 
in X . If not, the term 1 — becomes important and reflects the symmetry and periodicity of the 
lattice. 



The charges q x interact via a short-range (dipolar) potential 



U-{x-x') 



J_ V 



(2.12) 



N — > oo, \x — x'\ fixed, 



(2.13) 



where d(x,x') is the length of the shortest path between x and x', taking into account the periodic 
geometry. Furthermore, one obtains from the calculation that both potentials have the symmetry 
properties 

U±(X) = U±(X + N), 
U±(X) = U±(-X). 

Because of the range of the interactions, we will also call the long-range Coulomb interaction between 
the charges q x strong interaction and the short-range (dipolar) interaction between the charges q~ 
weak interaction. In the large N limit, for convenience, we rewrite the Coulomb interaction term in 
( 2.1C ), using ( 2.11 ) and charge neutrality, as 



N 



N 



x,x—l 



71-2 1xli'U+(x-x') = -7T 2 Y \ X ~ X 'kt<lx' - ~Jf [ J2 X <lt 

x,x' — l 



N 



(2.14) 



Inserting this expression in (2.1C) and writing out also the interaction potentials U p explicitly, we 
obtain eventually 



Tip — 



N N 

^(« + |Ec+E^ 

x—l x — l 

2 



r- N - 

^ 2 £ (3-2V2) d ^\-q-, 

x.x' — l 



Itt 2 N (m 



271 



47r m 



N N 



(2.15) 



E xq£ - 71-2 E \x-x , \q£qi>, 

x—l x,x' — l 



valid in the large N limit. The task will now be to minimise 7i P with respect to the four variables qg", 
q x , m, and n in order to find its ground state energy. This will be done in section III. 



5 



2.2 Antiperiodic boundary conditions 



Passing from PBC to antiperiodic boundary conditions (APBC) means changing the sign of the 
two horizontal bonds that belong to the plaquettes (N,l) and (N, 2). Under this change frustrated 
(unfrustrated) plaquettes remain frustrated (unfrustrated) . Thus the only modification needed to get 
the Hamiltonian 7i AP for APBC is to replace ft/N,i) by f(jv,x) + n m the first term in equation ( 2.15| ), 
i.e., to add h in the expression between parentheses in that term. 



2.3 Reflecting boundary conditions 

One obtains the Hamiltonian Tt R for the XY spin glass on an N x M lattice with reflecting boundary 
conditions (RBC) in the horizontal direction and PBC in the vertical direction by replacing the 
horizontal interactions in one single, but arbitrary column, say AT, by 

(<j)i + <j)j - 2-KTlij ~ TTij) 2 ■ (2-16) 

This amounts to reflecting the spins on one side of this column with respect to the reference axis. The 
ensuing modifications in passing to the Coulomb representation result in || 



{q r } n,m 



T 12 + 7T 21 
N + N 



mod 2,0 (2.17) 



with 

H R = vr 2 ^q r qr'U R (r-r') , (2.18) 

r,r' 

and r = (x, y) labels the centres of the plaquettes on the N x M lattice as before. We do not recall 

the explicit, general, form for the potential U R here, but rather use the variables and q~ defined 



in (2.9) and give the explicit expression of 7i Tt in the case of the tube lattice: 

jV JV 
H R = 7r 2 ]T q+q+,U+(x-x') + n 2 £ gjg", U~(x - x') (2.19) 

x.x' — l x.x' — l 

where 

U+(X) = f(l-^), JV-.oo.J^fixed.OgMgl, 

(2.20) 

TJ-{x-x') = ^(3-2V2) d ^^a(x,x') , N -> oo , \x - x'\ fixed , 

where a{x, x') = — 1 if the shortest path between x and x' crosses the bond irj^ (or tt 2 ^), and a{x, x') = 
1 otherwise. The interactions U R and U~ differ furthermore in their symmetry properties from those 
of PBC in that one has now 

U±(X) = -U±{X + N), (2.21) 
i.e., antiperiodicity of the interaction potentials. 
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2.4 The Hamiltonians in terms of electric field energy 



In this subsection, we rewrite the strong interaction part of the Hamiltonians (2.15) and (2.19) in 
terms of an electric field E T : 



(2.22) 



E x is the electric field between x and x + 1, and Eq is a constant background field whose value will 
be set later in such a way that the volume sum of the energy density E x gives the Coulomb energy 
of the Hamiltonians. The advantage of this rewriting is clearly seen in Appendix A when properties 
of the ground states of the Hamiltonians (for large N) are proven: E x is a local variable, whereas 

N 

E q x >Up R (x — x') involves all columns of the lattice. To determine the effects on the system's energy 

x' = l 

when a charge q Xo is changed in a given configuration is much easier in terms of the local variable 
E Xo , as is manifest in Appendix A. 



Squaring ( 2.22 ) and summing over x, we get 



jV N x N x x 

E El = NEq +2E q J: £ q+ + E E E lUt> 

x — l x — lx' — l x — lx' — lx" — l 

N N N 

= NEq +2E J2(N-x + l)q+ + £ £ [N + 1 - max^', x")]q+q+, 

x — l x' = l x"=l 



Using the identity 



max(x', x") 



\x'-x"\ x' + x" 



and rearranging terms leads to 
N 



N 

E^ 

x=l 



N N 

£o 2 + 2£oE<z, + + (E^ 2 

x—l x—l 

E E ^itit- ■ 



N 



N 



2E + £(x-l)g+ 



x' = l x" — 1 

The value of the constant background field Eo is obtained by setting 

JV N 

2n^E 2 x =- 2 E «Ut-U+ R (x - x>) 



(2.23) 



(2.24) 



(2.25) 



(2.26) 



In the large N limit, this amounts to comparing the expression in ( 2.25|) with the Coulomb interaction 



(2.27) 



part of Q2.15| ) and (|T|) [after insertion of ( |2.20D ]. This gives 

for H P and 7i AP , 



E 



2tt ■ 
N 

\ E it > for 1 

x=l 



and the Hamiltonians read hence 

JV _ _ 

H P . AP = 2ir 2 J2 E x+ IT 71 " 2 E (3-2V2) d ^ x kxQ x ' 

x=l 



N 



x^x' — l 
, N N 

x — l x—l 



2 "AP j 



N _ _ 

2^ 2 E^ + ^ 2 £ (3-2V2) d ^K(x,x')q-q-, 

x—l x,x' — l 



N 



(2.28) 
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where 




for PBC , 
for APBC , 



(2.29) 



with, from and ( fTf ) 



JV 

Qx = for Hp and H AP , 

x=l 



N 



it mod 2 = mod 2 for ft R 



(2.30) 



Having established the Hamiltonians for the different boundary conditions, we are now ready to 
determine those properties of their ground state configurations that are sufficient to calculate the 
typical energy difference between the ground state energies for N — > oo. 



3 The ground states for the different boundary conditions 

We summarise the problem to which the preceding sections have led. Each of the three expressions 



( 2.2S ) should now be minimised with respect to the variables n, {E x }, and {q x }. The {E x } are 



defined in terms of m and {q x } by (2.22) and (2.27), and the {q x } and {q x } are defined in terms of 



the original charges {q r } by (2.9). The variables m and n are integers, and the {q r } are half-integer 



or integer according to whether the plaquette r is or is not frustrated. 



The ground states of the Hamiltonians ( 2.2S ) possess in the large N limit the following properties: 



1. The electric field satisfies \E X \ ^ ^ for x = 1,. . . ,N. The constant background field takes values 

2. The charges q( x ,i)> Q(x,2) take the values 0,±^. 

3. On doubly frustrated columns, the charges <7( x ,i) an d Q(x,2) are equal if and only if \E x —i\ = |. 

These properties of the ground states, for N — > oo, are proved in Appendix A. Property 1 indi- 
cates that the system's ground state is within the set of states that minimise the electric field (i.e., 
the Coulomb) energy. For PBC and APBC charge neutrality implies En = Eq (cf. equation ( [2.22; )), 

N 

whereas for RBC antiperiodicity of the strong interaction potential leads to Ejy = —Eo = 5^5+ 

x = l 



(from equations (2.22) and (2.27)). Let us just note here that it was conjectured in |2| that, in the 
ground state of frustrated XY spin systems, the chiralities (charges) q r are likely to be zero on non- 
frustrated plaquettes and take the values ±4 on frustrated plaquettes. The above property 2 shows 
that, for large N, this is indeed the case for the tube lattice. 



We will now construct the set of states that have properties 1 and 2: We start placing charges 
1(x,y) — Oji^ successively from x — 1 to x = N, minimising the local electric field energy density 
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Figure 2: Construction of a state that possesses properties 1 and 2. We start placing charges Q( x ,y) from the left to the 
right. We proceed in a way that the local electric field stays as close to as possible. The signs indicate plaquettes 
with charges ± i . Plaquettes with no sign are charge-free (unfrustrated) . 



E x at each step and knowing that it changes by a half- integer amount (with q+ = ±i) on a column 
with one frustrated plaquette and by an integer amount (with q+ = for E x _\ = and qg" = 0,±1 
otherwise) on a column with both plaquettes frustrated (see figure 2). 

Having thus obtained a state that has properties 1 and 2, we see that it is always possible 
to partition the nonzero charges qt x ,y) into dipoles as in figure 2, grouping together two successive 
charges of opposite sign along the x-axis such that outside the dipoles the electric field is zero. [There 
is one exceptional case (see figure 3): for Tt K and Eg = i the last and the first charge placed are 
not part of any dipole, but of the same sign to take account of the antiperiodicity of the potential 
U£(X) = —U+(X + N), respectively En = —Eo, as mentioned above.] As one sees from figure 2, 
e.g. from the dipole containing charges on the columns x\ and x'i, dipole reversals do not change the 
Coulomb (i.e., electric field) energy of the system. So, the ground state of the system is found within 
a set consisting of chains of dipoles, satisfying properties 1 and 2, degenerate in Coulomb energy. The 
possibility of columns with two identical charges leads to the partition into dipoles not being unique. 
Property 3, however, introduces a further constraint on the set among which one finds the ground 
state. Furthermore, one can easily convince oneself that this latter property implies that one can 
reach every state that satisfies properties 1 to 3 from any other such state, by reversals of dipoles, 
for any given partition. In particular, the ground state of the system differs, for large N from a state 
as constructed above (with the additional constraint from property 3) by a reversal of dipoles. 

The exact ground state configuration remains unknown, but we know that it minimises the 
Coulomb energy independently of the other energies involved and have characterised the set of 
Coulomb energy ground states. [Let us just point out here that, when passing from PBC to RBC, 
one conserves by virtue of ( 2.27 ) and ( |2.3C| ) the property Eq — or E ^ for the ground states at 



both boundary conditions, due to the fact that ^ N 2 J N stays the same [see also Appendix A, espe- 



cially equation (A.9)]. So the Coulomb energy stays indeed the same.] The only remaining degrees 
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Figure 3: Example of a state with properties 1 and 2 for RBC, while Eq 



For RBC, proceeding in the construction 



(2.27), E 



as described in the text, the last charge to be placed is of the same sign as the first. This is obvious from equation 

JV 

\ ^2 which indicates that, in general, there is a surplus of two charges q( x ,y) with opposite sign to 

x=l 

Eq = ±5 (here Eq = i), to take account of the fact that for RBC the potentials Uj^ are antiperiodic (see (2.21)). But 

still all charges, but the first and the last one, can be grouped into dipoles as announced in the text. 



of freedom in this set are the directions of single dipoles: The degeneracy is lifted by the other terms 



in the Hamiltonians of equation (2. 28), which fix these directions. The amount of Coulomb energy 
of the system with the different boundary conditions being the same, it is the effect of these other 
terms that give rise to the difference between ground state energies when one varies the boundary 
conditions. We address this issue in detail in the following section. In spite of the fact that we ignore 
the exact ground state configurations, the above properties s! uffice to analyse and determine th e 
ground state energy differences for N — ► oo. 



4 Boundary conditions and ground state energy differences 
4.1 Generalities 

There is a general relation between the finite size scaling of the energy difference, AE^ ~°° JN^ 
(where J is the energy scale), of the ground states of a system for different boundary conditions and 
the corresponding correlation length, £(T), at a finite temperature T: The correlation length £(T) is 
set by AE& ~ k B T, hence 

[Let us just note here that the energy difference may be either concentrated in a domain wall or 
associated with a continuous variation of the order parameter.] 

In the Villain model, we may study the spin-spin correlation and the chirality-chirality correlation 
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by applying APBC and RBC §,§,§. So wc have to calculate 

where 4^, 4r\ and are the ground state energies under P, AP, and R boundary conditions, 

respectively. 

In this section and in Appendix B, we will call a dipole with charges Q( x ,y) an< 4 Q( x ',y') slanted, if 
y ^ y', and horizontal, if y = y' . [With this definition, the slanted ones include the vertical dipoles.] 
As the probability for a plaquette to be frustrated is the same for all plaquettes, a dipole is as likely 
to be slanted as horizontal. 



Furthermore, writing U~ for U~ at PBC ( [2.12D and for U~ at RBC ( pT20| ), the weak (dipolar) 
interaction U~(£) has the property 

oo 

U-(£)> U ~( £ ') ( 4 - 3 ) 

£'=e+i 

for all boundary conditions, so that we may approximate its effect by restricting the interactions of 
each nonzero charge q~ to those with its two nonvanishing neighbouring charges. Upon renumbering 
the nonzero charges q~ on the frustrated columns by a new index s — 1,2,..., N c (where N c is the 
total number of the frustrated columns with nonzero q~), we can finally rewrite the effective weak 
(dipolar) interaction as 

U~(x -x') = J2 U s q-q- +1 + U(0) ^(isf , (4-4) 

x,x' s—1 s—1 

where qN c +i = qi for PBC/APBC and qN c +i = —<h f° r RBC. The charges q~ take the values ±^ and 
±1, and the U s are independent quenched random interaction constants. Since in the set of states 
that we consider the Coulomb energy is boundary condition independent, we will deduce the energy 
differences, and AE^ (N — > oo), from (4.4) and from the global spin wave term in ( [2.2§| ). 

The large N limit is self understood in what follows. 



4.2 Antiperiodic boundary conditions 



The actual ground state minimises the second and third term in TC r and Ti Ar , equation (2.28), within 
the space of degenerate ground states of the Coulomb energy, characterised in the preceding section. 
The second term, rewritten in equation (4.4), is the weak interaction between the q~ . In Appendix B 
we show that its lowest-lying excitation lies typically an energy amount ~ JN~ Vc above the ground 
state, and determine the exponent y c — logs (3 + 2y2) = 1.7972 .... The third term is 



in Hp , 



(4.5) 
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where we have used the neutrality condition to write J^Qx — 2 X Q(x,x)- The terms in ( |4.5| ), of order 

X X 

iV -1 , have their origin in a global spin wave, i.e., of wave length > N, which helps the system to 
adjust to its boundary conditions when there is a rotational mismatch (cf. [|[ ||). 

One might wonder if it is always possible, by choosing n properly, that the terms in ( |4.5| ) vanish 

N 

in the ground state in H P and/or H AP . As n + G {0, ±A,±1,...}, this depends obviously 

ce=l ^ 

on the number of nonzero charges ?(x,iV Hence we have two cases: 

(i) the number of frustrated plaquettes (x, 1) is even; 
(ii) the number of frustrated plaquettes (x, 1) is odd. 

Correspondingly: 

N 

\ E Q(x,i) e {0,±i,±l, . . .} case (%) , 

(4-6) 

^?M£{44--} case (nj . 



!E=1 

We investigate these cases further. 



(i) Even number of frustrated plaquettes (x, 1): Given the set of ft( x ,i) and possibly reversing a sequence 
of dipoles as in Appendix B to get | ^ 9(x,i) + E ~~2^" integer for PBC and halfmteger for APBC 

X X 

or vice versa, the terms in ( |4.5| ) vanish for both Hamiltonians by a proper choice of n. As there is 
a difference of \ in the term in parentheses in (4.5), the ground states of Tt P and TC AP differ by a 



reversal of a sequence of dipoles containing an odd number of qt Xt i) ^ 0, i.e., a sequence of dipoles 
among which an odd number is slanted. Hence we obtain, reinserting J, 

A£^ip ) - ±JN~ V % case (i) , (4.7) 

where the sign indicates that either state, at PBC or APBC, has the lower ground state energy. 



(ii) Odd number of frustrated plaquettes (x, 1): Here, from (4.6), the terms (4.5) in TL P and Tl AP are 

2 

always nonzero and the optimal value of n will give an energy Jfjy irrespective of the directions of 
the dipoles. These will thus be determined by the weak interaction only and be the same for both 
boundary conditions. The energy difference is hence 

AE AP ] = 0, case (ii) , (4.8) 

in this case. 

4.3 Reflecting boundary conditions 

The ground states at PBC and RBC minimise the second and the third term in TL P and the second 



term in 7i R , equation ( 2.28 ), within the set of states characterised in the preceding section. So, again, 
one has to distinguish between an even and odd number of qi x ,x) 7^ 0, when calculating the typical 
ground state energy difference. 
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(i) Even number of frustrated plaquettes (x, 1): We saw in section 4.2 that the term for PBC in fl4.q ) 
vanishes in the ground state of Hp and that, for half of the samples, one finds the ground state at PBC 
by minimising the weak (dipolar) interaction. If one changes to RBC for this half of the samples, while 
keeping the PBC ground state configuration, one will almost always be able to lower the energy [Note 
that the sign of the bond that passes column N, which is almost never the weakest one, has changed 
(!)]: one just has to reverse a sequence of dipoles starting at column N in such a way that only one 
of the weakest bonds is broken. For the other half of the samples, the cancellation of the global spin 
wave term implies a reversal of a sequence of dipoles in the configuration after minimisation of the 
dipolar interaction. For these latter samples, when one changes to RBC, one has to reverse again the 
same sequence that was reversed to get the ground state at PBC. In both cases, this leads to 

AE { n N) ~ ±JN- y % case (i) , (4.9) 

for the difference in ground state energies. The minus sign applies for the first half of the samples, 
the plus for the second half. 

(ii) Odd number of frustrated plaquettes (x, 1): The global spin wave term never vanishes in PBC, so 
that 

2 

AE { R N) ~ ^JN- 1 , case (ii) , (4.10) 

neglecting a possible contribution of order JN~ Vc . 



5 Conclusion 

We have studied the XY spin glass with ±J bonds on a tube lattice. This system has both a 
continuous (spin) and a discrete (chiral) symmetry, and hence two order parameters play a role. Our 
purpose was to determine the divergence, for T — > 0, of the chiral and the spin correlation lengths, 
via the finite size scaling of the ground state energy differences under different boundary conditions. 
In the presence of two symmetries, the usual single-symmetry relation between the finite size scaling 
exponents of the ground state energy difference and the correlation length has to be extended in a 
nontrivial way. Nevertheless, the spin correlation length exponent y s (see equation ( |4.l| )) is given 
by the energy difference when one passes from periodic to antiperiodic boundary conditions, namely 



(AE^)) 2 ~ JN Va . New boundary conditions, reflecting ones, were introduced || to determine the 
chirality correlation length exponent y c . 

The difficulty in performing such an analysis on a general N x M lattice is that one does not know 
how to construct the ground states of the disordered systems. The tube lattice, of this work, however, 
just as the ladder lattice studied earlier allows for a precise theoretical analysis of this relation. In 
contrast to the ladder lattice, the tube lattice still has long-range interactions between its chiralities, 
and is therefore closer to a two-dimensional system. 

We first apply the well-known transformation |9) of the XY spin glass into a Coulomb gas, a 
system of chiral variables (also called charges) . The resulting effective Hamiltonian can be cast in the 



form (2.28) where it is the sum of three terms: 
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(i) A one- dimensional Coulomb interaction, linearly increasing with distance, between charges 



q>2 , . . . ,qjf; in ( 2.26 ), this term has been expressed as the volume sum of the energy density of 



the electric field E x . 

(ii) A "dipolar" interaction that decreases exponentially with distance between the q± , . . . , qjf. 

(Hi) The energy of a spin wave needed to match PBC or APBC (but absent under RBC), and whose 
wavenumber depends on the total electric dipole moment. 

The third term disappears in the thermodynamic limit. Its relevance for a finite size scaling 
analysis was first pointed out by Fisher, Tokuyasu and Young JlO[ ] . Moreover, the three terms are, on 
the one hand, coupled by local constraints, that link the allowed values of q£ and q~ with the fixed 
values of the ferromagnetic or antiferromagnetic bonds 7Ty between the spins on the lattice, and, on 
the other hand, by a global constraint on the total charge (zero for PBC and APBC, and even or 
uneven for RBC). Taking these constraints into account, we identify the low-lying excitations of the 
three terms, respectively: 

(i) Coulomb excitations that cost an energy of order J. 

(ii) Chiral excitations, obtained by reversing a sequence of chiral variables, that cost an energy ~ 
JN~ Vc withy c = 1.7972... . 

(Hi) Global spin waves that cost an energy ~ JN~ X . 

The ±J XY spin glass on the ladder lattice || consists of both interactions (ii) and (Hi). Due to 
the additional long-range interaction (i), the tube is closer to the two-dimensional model. 

In spite of the number of interactions in competition, we were able to characterise and delimit the 
set of charge configurations, within which lies the ground state. In the configurations contained in 
this set, the charges take the values ±i on the frustrated plaquettes and zero on the others, and form 
a chain of dipoles. 

We now give a summary of our results, and recall numerical results for comparison. When changing 
boundary conditions from PBC to APBC, or RBC, it is the excitations (ii) and (Hi) that give both 
energy differences, AE^' and AJ?r . This implies the same conclusions as in ||: Firstly, the 
ground state obtained with P boundary conditions can adjust to AP boundary conditions via a chiral 
excitation, so that 



(AbIp') 2 '-^, y c = 1.7972 ... . (5.1) 

The last equation contains no reference to spin waves and means that y s = y c . Secondly, passing 
from P to R boundary conditions releases a global spin wave (as was first observed by Kawamura and 
Tanemura || in d = 2) in half of the samples, but does not do so in the other half. 

In d = 2, Kawamura and Tanemura performed a numerical analysis of the different ground state 
energies of the cosine XY model. They find, as N — > oo, 



(a4p') 22 =s aN~y° , y s » 0.84, 

j. (5.2) 

( e W_ e W)2 2 ~ bN-y + 0(N-*), y c « 0.38, 
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where a and b are constants, and a new quantity, namely 



(5.3) 



has been introduced. Thus, they get two distinct exponents y s and y c , with y s > y c , and conclude 
that the chiralities order on a longer scale than the spin variables. 



For the tube, upon collecting our results [equations (4.7),(fD]),(4.9) and ( 4.10p ], we get 



(4 J 



(N) (N) 



) 2 = — JN- 1 + (D(N- y c) , 
8 



(5.4) 



i.e., the R boundary conditions probe a global spin wave term proportional to TV l . If we now 
conjecture on the extrapolation of our results to d = 2, then we expect for the quantities of equation 



(5.2) that AE Ar would yield a chiral exponent y c as in (5.1) but with a smaller value (since y c should 
vanish at some, still higher, lower critical dimension); and that e R would yield the spin wave exponent 
d — 2 = 0. Instead, in contrast, Kawamura and Tanemura interpret their simulation according to 
(|5.2|). We expect that simulations on larger 2d systems will confirm our scenario. 



Appendix A 

In this appendix, we prove that in the large N limit the ground states of the system for the different 
boundary conditions possess the properties 1 to 3 announced in section III. In the calculations, we 
write the expression of the weak (dipolar) interaction in its form at PBC/APBC (equation ( 2.12p ), 



with again J = 2. The arguments are nevertheless readily rewritten for RBC, including the appropri- 
ate factors of a(x, x') (equation ( |2.2C| )). Furthermore, we neglect the global spin wave term O(j^) that 
appears for PBC/APBC. Upon proper choice of n, this term contributes at most ™ to the ground 
state energies at PBC/APBC, which is small for N — > oo, in comparison to the other energies involved. 



In preparation of the proofs of the ground state properties 1 to 3, we show, in a first step, that 
(i) in the ground state, the charges q~ take values \q~\ ^ |, 
and, using (i), in a second step, that 
(ii) in the ground state, the charges q~ take values \q~\ S 1. 
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Proofs of (i) and (ii): 



(i) In the ground state, the charges q x take values such that \q x | ^ 1. 
Let us look at some state with charges q~' such that 

<Z = max{|^°|}>2 . 



(A.l) 



Let q x f be a charge with \q x f\ = 
positive without loss of generality. 



q. For reasons of charge reversal symmetry, we may take q x 



Consider now the state in which q~^ is changed into q~^° — 2, while all other charges q~'° and all 
charges q x '° are kept unchanged. [Note that, by an appropriate change of q< Xt x) an d Q(x,2)i one can 
add an arbitrary multiple of 2 to some charge q~ without changing g+ (see ( |2.9[ )).] The difference in 
energy AE between the (final) state, with q~>° changed, and the initial state can readily be calculated. 
As the charges q x are unchanged, it comes from the difference in weak interaction energy only. For 
q x f ^ 2, one finds 



AE 



f- 2 [(«■ 



2) 2 -q 2 }+^ 2 [«-(«- 2)] 2 £(4 



s>0 



x-\-s 



q-l° s )(3-2V2Y 



-i(q - 1) + 4 £ (9-4.° + £1° )(3 - 2^2)' 



s>0 



(A.2) 



We have 



4^(C4-°+Cl°)(3-2V2>' 

s>0 



E 

s>0 



(3 - 2\/2>' 



and thus, summing the geometric series and using q ^ 2, 

n/2 



-4(g-l)+4(V2-l)( 



< . 



(A.3) 



(A.4) 



So the final state is lower in energy than the initial state. Hence a state, in which q ^ 2, is not the 
ground state. 

(ii) In the ground state, the charges q~ take values such that \q~\ ^ 1. 

Let us take some state with charges q x '° — 0, ±i, ±1,±|. Be n the number of charges ±|, with 
n > 0. Suppose for the moment that n < N. There exists hence a sequence of charges q X( : , q. 



210 + 1'' 



if (1 ^ 
1^1). 



/ < 



that is enclosed between charges q x ,0 of absolute value ^ 1 (i.e., 



-,o 

"xo+n'-l ~ 

I I I _ >° 

\Qxo-l\>\"x +n 

Consider now the state in which all charges q x '° — ±| in this sequence are replaced by =pi, while 
all other charges q~'° and all charges are kept unchanged. As in (i), the difference in energy 
AE between the (final) state, with the changes, and the initial state is due to the difference in weak 
(dipolar) interaction energy only. The difference AEq in weak interaction energy, coming from the 
self-interaction terms in the sum, is 



AE = -2^p7rV , 
the one from the nearest-neighbour interaction terms is of absolute value 



A^j <2^tt 2 (3-2V2) 



("'-1)7 + 



2 .|, 



2^ 2 (3-2V2) 



(n'-l) 7 + 



2 -2 a 



(A.5) 



(A.6) 
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while the energy difference of all other terms is of absolute value 

|Ai^4^ £-2^9 ^ (3-2V2)' 3 n> 

1 8 l-(3-2V5)4 8 l-(3-2V2)4 1 y 

To obtain the last inequality, we have substituted for all charges q~ , but the ones in the sequence, 
the maximal possible absolute value | and taken all terms in the sum to be negative in the initial 
state and positive in the final state [which is obviously an upper bound for the energy change, but 
impossible to realise] . The overall energy difference is thus 



AE £ -2^ttV + ^tt 2 (3 - 2V2) [bin' - 1) + 8] + 12^ 2 ^f^ n' 
= -^ttV(29 - 20V2)-9 + 6V2] <0 . 



(A.. 



Again, the energy of the final state is lower than the energy of the initial state. If n = iV, a similar 
reasoning (with n' = n = N) leads to the same conclusion. This proves (ii). 



We are now prepared to show that the ground state has the properties 1 to 3, announced 
in section III. 



1. In the ground state, the electric field satisfies \E X \ = |. The constant background field takes values 
We first note that the electric field is constrained to be of absolute value ^ 1 for the ground state 



configuration. This can be seen as follows. From ( [2.28D , we see that minimisation of the Coulomb 
interaction part of the Hamiltonians means minimising the mean square value of the local electric field 
E x . This leads for all Hamiltonians to 

if *n+*n mod 2 = 0, 
E ={ ' (A.9) 
±\, if N ^ N mod 2 = 1, 

for the constant background field Eq. The electric field, E x , should be locally optimal, that is stay 
as close to as possible. For any given state with given sets of charges {q x } and {q x }, one can go 
successively through the system, from x — 1 to x — N, changing, whenever necessary, the charges q x 
in such a way that the value of the electric field is bounded in absolute value by 1 in the final state: 
At every column with no frustrated or two frustrated plaquettes, the change in electric field can be 
bounded to be 0,±1 by an appropriate change of g+, if necessary, and to be ±| on every column 
with exactly one frustrated plaquette. [By the definition of q x and q~ , equation ( |2.9| ), one can add an 
arbitrary multiple of 2 to some charge q+ without changing q~ , by an appropriate change of q(x,i) an d 
9(x,2)-] Whenever one encounters a column xq, during the above procedure, where the electric field 
jumps to a value \E Xo \ ^ |, one changes q XQ by the appropriate amount, as well as the next nonzero 
charge qj, say at x' , by the opposite amount (to conserve charge neutrality). By the definition of the 
electric field, cf. equation ( 2.22| ), every time one changes the charges at xo and x' only, the electric 
field remains unchanged on columns x < xq and x ^ Xq. So in the end, while keeping the set of 
charges {q~} fixed, the value of the electric field is bounded in absolute value by 1. In particular, the 
final state is lower in energy than the initial state. 
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Let us now take a charge configuration of the system (with charges q x ,Q , <j£~' ) such that the electric 
field jumps, say at column xq, to a value l-E^I = 1, and stays at this value until column x' (> xq), 
where one finds hence the next nonzero charge g^ ■ Consider the state in which q£'° is changed by an 
amount of 1 to some value \E'\ ^ | and q£,' by the opposite amount (to conserve charge neutrality), 
while keeping the charges qp° (and q x '°) on all other columns fixed. We observe that, from ( [2.22j ), 
the electric field is unchanged for x < xq and x ^ x' and that by the definition of q£ and q~ , the 
charges q~>° and q~,' changed. Anyhow, from (ii), the absolute value of the charges q~ can be taken 
to be bounded by 1 in both states, as otherwise the energy of the state can still be lowered. The 
difference AE F in strong interaction energy between the final state and the initial state is 



AE l 



-2tt 2 - \x q -x' \ , 



(A.10) 



while the difference AE* in weak interaction energy is bounded by 



|A£^| 



< 



2^7T 2 

o 



s>0 



1 



s>0 c 



,0 
+s 



.1 .-, 



J(3-2V2>' 



(A.ll) 



[The factors 2 in the last inequality stem from the fact that both the energy of the initial state and 
the one of the final state enter in the difference.] Hence we get for the total energy difference 



A£<-2tt 2 -+4— tt 2 [1 + 2(V2-1)] < -^^—V , (A.12) 
4 8 2 

1. e., the energy of the final state is lower than the one of the initial state. Thus \E X \ ^ ^ for the 
ground state configuration. 

2. In the ground state, the charges <7( x ,i)> Q(x,2) take the values 0, ±i. 

From (ii) and 1 , the absolute values of the charges on the plaquettes in the ground state are bounded 
by 1. In the case that E x _i = in the ground state, it is easy to see from (ii) and 1 that the charges 
q<x,i): q(x.2) on the plaquettes of column x are of value 0,±i. 

Let us thus consider the case = \. If there is at least one frustrated plaquette on column 

x, it is again obvious from (ii) and 1 that q( x ,i) an( A 9(x,2) are °f absolute value ^ h. Suppose now 
that there are two nonfrustrated plaquettes on a certain column xq . From what we have stated at the 
beginning of this paragraph, it could be that there is a charge ±1 on one of the plaquettes. Let us 
furthermore suppose that there is a sequence of columns xq , Xq + 1,- • •, xq + n — 1 that carry charges 
q~ of absolute value 1 . 

If the number n of charges q~ in the sequence is even, consider the state in which all the charges 
Qx >Qxo+i> ' ' 'Axn+n-i are changed to [this is possible as it conserves charge neutrality, as one 
convinces oneself with the help of 1 ] . The absolute value of the electric field is the same in both the 
initial and the final state, so that the difference in energy is again due to the weak interaction only: 



. ~ v2 9 d V2 9 3 — 2v2 

AE < — - — 7i Ti + 4—ir 2 n — 

~ 8 8 l-(3-2V2) 



^7r 2 n[2V2-3] < 
8 



(A.13) 



and the energy of the final state is lower than the one of the initial state. 
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In the case that n is odd, one can take n = 1 and l^-il'kx+il = \ without loss of generality, in 
view of the preceding paragraph. Suppose first that one of the charges q~_ 1 ,q~ +1 is of absolute value \ 
in the ground state, say and say [(ftx+inl = \. Let us compare the energy of this state with the 
one in which q~ is changed to and q( x +i,i) replaced by — q( x +i.\) [note that q~ +q( x +i,i) = — 9(x+i,i) 
from 1, so that the absolute value of the electric field remains unchanged and charge neutrality is 
conserved]. Using (ii), the difference in energy is 

A£ ^ f tt 2 [-1 + 2(\q-f iq -'°\ + IC'V+il + l^+k~+ 2 l)( 3 - 2V^) 

+ 2\ - a°ia°2i(3 - 2V2) + 2 e iffx ,0 ikx-°. + <a° i( 3 - 2 V2) s 

S>1 

+ 2 • 2 E la°illC4°i- s + <a°i +s K 3 - 2V2)^ 

^tt 2 (-17+ 12V2) < . 



S>1 



(A.14) 



-l + 8(3-2V2)^+8 (3 ~ 2V ^ )2 



2^2-2 



Again the final state is lower in energy than the initial state. Secondly, if g r ~_ 1 = q~ +1 = 0, there 
are again two possibilities. Either there are only nonfrustrated plaquettes on the columns x — 1 
and x + 1, or, on at least one of them, both plaquettes are frustrated, say at x + 1. In the last 
case, one has <?+ = -2g (a . +ljl) = -2g (x+li2) from i; hence changing q~ to and g( x+ i,i), ?(x+i,2) to 
~ — 1(x+i,2) conserves charge neutrality, the absolute value of the electric field and does not 
change q x _i = q x+1 = 0. Using (ii), the energy difference between the final and the initial state is 
then 

-1 + 2 E |<£ ,0 ll<fe-° + <a°l(3 - 2V2)* 

^tt 2 (-15 + lOy/2) < . 



AE < ^ 



< V2 



\/2^2 



4 (3-2v^ 
2/2-2 



(A.15) 



In the case that q( x -i.i) — q( x -i,2) = <Z(rr+i,i) = 5(x+i,2) = 0, there is a column x — s or x + s 
such that q( x - a ,i) — Q(x-a-.2) — 9(:e+<t,i) = q( x +a.2) = for all 1 < a < s and that one of the 
charges q(x-s.i)j < l(x-s.2)j < }(x+s,i)7 c l(x+s.2) is nonzero. Considerations analoguous to the ones earlier in 
this paragraph lead again to the conclusion that a state in which there is a nonzero charge on a column 
with two nonfrustrated plaquettes is not the ground state. 

So in the ground state, the charges (?( Xl i), q( x ,2) take the values 0,±i only. 



3. In the ground state, the charges q( x ,i) and <?( x ,2) are equal on doubly frustrated columns, if and only 

if\E x - X \ = \. 

Following exactly the same lines as in 2, one shows that a state in which q( Xo ,i) = —q( X() .2)i i- e -, 
\q~ \ = 1, on some doubly frustrated column x , while |-E;e _i| = \, is higher in energy than the 
one in which q( Xo ,i) — q( X[l ,2), i-e., \q Xo — 0| [the appropriate change of some other charge, as in 2, 
to conserve charge neutrality, is tacitly understood]. It is evident from 1 that, in the ground state, 
9(a;o,i) = — 1(x .2) on some doubly frustrated column xq, if E Xo -\ = 0. 
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Appendix B 



We calculate the typical energy change in a tube of length N, in the limit N — > oo, when a sequence 
of dipoles is reversed to adjust to APBC or RBC from PBC, as in section IV. This amount of energy 
is related to the typical length of the longest interval that contains no nonzero charges q~ between 
two such (nonzero) charges. The reversals that are effectuated in section IV do not necessarily involve 
the longest of these intervals because of the constraints: one must not brake up a dipole, or in certain 
cases has to reverse not just any sequence, but one that contains an odd number of slanted dipoles. 
Nevertheless, these reversals will still typically involve intervals that are of the same order as the 
longest interval. 

Since each plaquette is frustrated with probability i independently of the others, a tube of length 



N will typically contain -j- nonfrustrated columns, with one frustrated plaquette, and ^ where 
both plaquettes are frustrated. From property 3 of the ground state configuration, half of the doubly 
frustrated columns will typically contain charges that belong to the same dipole and the other half 
charges that belong to two different dipoles. This is due to the fact that at a given column E x = or 
E x = ±5 (i.e., q£ integer or half-integer) with equal probability; in the first case, one has q~ = ±1, 
and in the other, q~ = 0. So, there are more charges q~ = than there are nonfrustrated columns. 
Since it is typically half of the doubly frustrated columns that give q~ = 0, we sec that, again typically, 
|iV of the columns carry a nonzero charge q~ , while |JV of the columns are neutral in q~ . So the 
number of intervals between two nonzero charges q~ , that contain no other such (nonzero) charge, is 
|iV. The probability p(£) for the two subsequent nonvanishing charges q~ to be at a distance i (i.e., 
to be separated by an interval of I — 1 columns containing no nonzero charge q~) is 



p(£) 



= 1,2, 



(B.l) 



Let Pjv(m) be the probability distribution for the length m of the longest one of these distances. 
Obviously, P/v(to) equals the probability that all |iV intervals have length i ^ to, minus the proba- 
bility that they all have length I ^ m — 1 ; explicitly 



Pv(m) 



iN 



5>w 



(B.2) 



When TV is large, Pjv(m) will be peaked around some large value of to. Its scaling form can be 
obtained if one transforms from to to to' according to 



to = 7 log N + m , 
with 7 to be determined. Indeed, upon using ( |B.3| ) and ( |B . 2| ) , one finds 

with x = e~3vl) N 
This shows that P/v(to) is effectively nonzero only for argument values 

to = logs N + 0(1) , 



Pjv(to) N ~ CC x(1 — IE?), 



(B.3) 



(B.4) 



(B.5) 
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and that the appropriate scaling limit reads 



N -►oo, 

m! finite, fixed, (B.6) 

and 7 = t^t ■ 



Having thus obtained the typical length 

log TV 

m=-^ x (B.7) 
logf 

of the longest interval, we are able to determine the typical energy change of the mentioned dipolc 
reversals. From the effective weak interaction (4.4) we deduce that the typical energy change due to 
such a reversal is of order c U (m) , where 

U(m) = (3 - 2V2) m (B.8) 

is the energy change when breaking up a bond at distance m and the constant c is of order unity The 
exponent y c is then obtained from the defining equation 

cU{Q)N- y " = cU{m) (B.9) 

which gives 

y c = logs (3 + 2\/2) = 1.7972... (B.10) 
for the exponent of the chirality-chirality correlation length. 
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